Multi- Welled Tunneling Model for the Magnetic-Field Effect in Ultracold Glasses 
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Puzzling observations of both thermal and dielectric responses in multi-silicate glasses at low 
temperatures T to static magnetic fields B have been reported in the last decade and call for an 
extension of the standard two-level systems tunneling model. An explanation is proposed, capable 
of capturing at the same time the T- and _B-dependence of the specific heat C p and of the dielectric 
constant e in these glasses. This theory points to the existence of anomalous multi- welled tunneling 
systems in the glasses - alongside the standard two-level systems - and indications are given for 
glasses which should achieve larger electric magnetocapacitive enhancements. 

PACS numbers: 61.43.Fs, 77.22.Ch, 77.22.-d, 65.60.+a 
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The last decade has seen much renewed interest in the 
physics of cold non-metallic glasses, materials displaying 
some universal physical properties attributed to the low- 
energy excitations characterising most amorphous solids. 
The two-level systems (2LS) tunneling model (TM) Q 
has been rather successful in explaining a variety of in- 
teresting phenomena dealing with the thermal, dielectric 
and acoustic properties of structural glasses at tempera- 
tures T < 1 K. Limitations and failures of the 2LS TM 
(treating cooperative motion in terms of single particles, 
mostly) on the other hand, have been also discussed 

These materials do not present, normally, any remark- 
able magnetic-field response phenomena other than pos- 
sibly a weak contribution from trace paramagnetic Fe- 
impuritics. It came thus as a great surprise when mea- 
surements showed [3J that in some multi-component sil- 
icate glasses (but not in pure a-SiCh) one is able to ob- 
serve changes in the dielectric constant e(T, B) and al- 
ready in magnetic fields B as weak as a few Oe. A typi- 
cal glass showing a strong response has some 100 ppm 
Fe 3+ and the composition Ai2 03-BaO-Si02 (in short 
AlBaSi-O, in this paper). Measurements made on a 
thick sol-gel fabricated AlBaSi-0 film showed changes in 
Se/e = [e(B) — e(0)]/e(0) of order 10~ 4 and characterised 
by an enhancement peaking around 0.03 T for 10 mK 
< T < 200 mK then followed by a reduction of e for B > 
0.1 T. A further enhancement was also observed at much 
higher fields (B > 10 T). 

Another, cleaner, multi-component silicate glass 
(borosilicate BK7, with 6 ppm Fe 3+ ) and a dirtier one 
(borosilicate Duran, with 120 ppm Fe 3+ ) have shown 
similar - but weaker - magnetic anomalies, \Se/e\ ~ 
10~ 5 , seemingly excluding the paramagnetic impurities 
as their source Yet another multi-silicate glass, a- 
Si02+ K Cj,H 2 was investigated confirming the unusual 
findings in AlBaSi-O. A convincing explanation for the 
unusual electric magnetocapacitance behavior of these 
cold glasses has not yet been found. 

Recently, some consensus has been gained by the idea 
of a coupling of the standard 2LS to the magnetic field 



via nuclei in the glasses carrying an electric quadrupole 
moment as well as a magnetic dipole one Jf|. The nu- 
clear mechanism is supported by some features of the 
observed magnetic-field dependence of the polarization 
echo (PE) experiments in the mentioned multi-silicate 
glasses in the millikelvin range 0]. Moreover the ampli- 
tude of the PE in Glycerol glass was shown to become 
strongly B-dependent only upon deuteration and thus 
the introduction of quadrupole-moment carrying nuclei 
in the glass Q. However, though pure a-SiC>2 (devoid 
of quadrupole-moment carrying nuclei) shows no PE- 
amplitude magnetic field dependence Q , Glycerol glass - 
deuterated or not - shows no measurable ^-dependence 
in its dielectric constant . The nuclear approach is also 
unable to account for the magnitude and features of the 
B- and T-dependence of Se/e for the multi-silicate glasses 
10] and the (also unusual) B- and T-dependence of the 
heat capacity C p of the latter - not entirely linked to their 
Fe impurity contents [HI, HI] _ is not even addressed. To 
add to the mystery, the acoustic response (also linked 
to the 2LS coupling to phonons) of borosilicate glasses 
BK7 and AF45 has been found to be independent of B 
fl3l ] . Table I summarises this rather puzzling experimen- 
tal situation. Therefore, either the nuclear explanation 
is specific to the PE magnetic effect, or an entirely new 
explanation should be found for all of the observations. 

The purpose of this Letter is to begin to give a rationale 
to the situation in Table I, as well as to stimulate further 
experimental and theoretical research. A novel explana- 
tion, already shown to account for the unusual behavior 
of C P (T,B) in the multi-silicates fjjj], is here shown to 
explain the behavior of the dielectric constant as well. 
This simple theory is centered on the observation, in 
computer simulations and experiments [3], that multi- 
silicate glasses present in their atomic structure both the 
connected network of Si04-tetrahedra and a collection of 
pockets and channels of non-networking ions showing a 
tendency to form microaggregates and to partially de- 
stroy the Si04-network. Fig. 1 shows as an example a 
snapshot of a simulation of the (Na20)-3(Si02) glass il- 
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glass type 


Ref. 


SC P 5e 5v s 8 Ape 


a-Si0 2 




NO NO ? NO 


a-Si02+iC y H z 


[5] 


? YES ? ? 


AlBaSi-0 




YES YES ? YES 


Duran 




YES weak ? YES 


BK7 




weak weak NO YES 


AF45 


[13] 


? ? NO ? 


Glycerol 


0]J8] 


? NO ? weak 


d-Glycerol 


f9]J8] 


? NO ? YES 



TABLE I: Presence of magnetic-field induced variations in the 
physical properties of some cold glasses. C p : heat capacity, e: 
dielectric constant, v s : sound velocity, Ape'- PE amplitude 
[?: no investigation known] 



lustrating such situation. The present theory proposes 
that the magnetic effects arise from anomalous tunneling 
systems (ATS) forming in the cooling of such structure 
within the non-networking (or network-modifying, NM) 
pockets and channels, whilst the Si04-network remains 
the nest of the ordinary non-magnetic 2LS. 



and are taken, normally, to be distributed in the glass 
so that A and In Ao (roughly, the potential barrier) have 
a uniform distribution: V-2Ls(A,Aq) = P/Aq, P being 
a material-dependent constant. This description holds 
for the network-forming (NF) TS. For the NM ATS in- 
stead, the simplest 3D generalization of the 2LS TM is 
that of other fictitious charged particles, each moving in 
a multi-welled 3D potential (see Fig. 2) and coupling to 
the magnetic field through their orbital motion. For the 
simplest case of n w = 3 potential wells, one can use: 
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where Do — TiVte~ UB / hn is some 3D barrier's trans- 
parency (Ub barrier height and Q single- well frequency), 
Ei,E2, E3 are the single wells' ground-state energy asym- 
metries (Ei + E 2 + E 3 = 0, say) and 
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FIG. 1: Molecular dynamics snapshot of the structure of 
sodium trisilicate at 2100 K at the density p = 2.2 g cm" 3 . 
The blue spheres that are connected to each other represent 
the (NM) Na atoms. The Si-0 network is drawn by (NF) 
yellow (Si) and red (O) spheres that are connected to each 
other by covalent bonds (from [1J|, by permission). 

In order to couple the ATS to the magnetic field, a 
simple 3D generalization of the 2LS TM is in order. As 
is known [l| , in the TM the cold glass is thought to have 
few remaining degrees of freedom rappresented by fic- 
titious "particles" , each moving quantum-mechanically 
within a ID double-welled potential. At low temper- 
atures only the ground states of the individual wells 
\i) (i=l, 2) are relevant and in this representation the 
Hamiltonian of a single 2LS reads Ho — — \Ao z — iAo<7 x 
(ct m Pauli matrices) with A the ground-state energy 
asymmetry between the two wells and Ao the barrier's 
transparency. These two parameters are linked to the 
real potential's details as discussed in the TM literature 



is some Aharonov-Bohm phase for a "particle" carrying 
charge q tracing a closed path of area Sa threaded by a 
magnetic flux <&(B) ($ = h/\q\ = <po\e/q\ being the ap- 
propriate flux quantum (tpo being the elementary one)). 




(b) 





V(x,y) 















(c) 
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FIG. 2: 2D representation of the plausible source of magnetic- 
field sensitive (anomalous) tunneling systems in (e.g.) the 
AlBaSi-O glass. The tight vitreous-Si02 structure (a) is bro- 
ken up by the Al- and (large) Ba-atoms (b), thus leaving many 
metal ions free to move as "particles" in a 3D tunneling po- 
tential characterised by n„-minima, with n w > 3 (c). 

In this model the choice Do > can be thought of 
as arising from the coherent tunneling motion of a small 
cluster of NM-ions; this will lead, as is seen, to high values 
of the product \q\S A D . For D ee y/Ef + E\ + E% < 
Do and weak fields ((j> <C 1) the lowest energy gap 
of the model is seen to open with increasing magnetic 
field according to the simplified expression 12j AS ~ 
y/ Dq(/) 2 + D 2 containing the main physics of this model 
(regardless of the value of n w ). One more assumption of 
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this theory (seen to explain the nearly-flat T-dependcncc 
of C p for B — in some temperature range for these 
glasses [lH) is to take the parameters' distribution uni- 
form for Ub, but favoring near-degeneracy (to a degree 
fixed by a lower bound D m i n ^ for D) for the {Ei}: 
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This anomalous distribution for the ATS can be thought 
of as arising from a degree of devitrification in the mate- 
rial, as measured by the parameter P*. In fact, it is 
reported 15j that thick glass films prepared with the 
sol-gel technique (such as the AlBaSi-0 films of the ex- 
periments) are multiphase materials with microcrystals 
embedded within an amorphous glassy matrix [161 1 1 71 ] - 
Indeed, amorphous solids with the general composition 
Al2 03-MgO-CaO-Si02 are termed "glass ceramics" in 
the literature [l6| owing to partial devitrification occur- 
ring. It seems thus reasonable to imagine that in these 
glasses some NM-ions provide nucleation centres for the 
microcrystals and that, therefore, TS presenting near- 
symmetric wells will be favoured in the neighborhood of 
and inside such crystallites. 

This approach has provided a good description [l2| 
of the C P (T,B) data for AlBaSi-0 and Duran [Hi; one 
can treat the ATS as effective 2LS having gap AS for 
"weak" fields. Within this picture, the linear-response 
quasi-static resonant contribution to the polarizability is 
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contains the single- well dipoles Pi = qa^. This expres- 
sion assumes vanishing electric fields and no TS-TS in- 
teractions, a situation which does not wholly apply to 
the experiments. To keep the theory simple one can still 
use Eq. ((4]) and the analogous one for the relaxational 
contribution to the polarizability. Eq. ^ must be aver- 
aged over the random energies' distribution Q ([. . ] av , 
responsible for the high sensitivity to weak fields) and 
over the dipoles' orientations and strengths ((. . .)). For 
a collection of ATS with n w > 2 this averaging presents 
serious difficulties and one must resort to the decoupling: 



G^8{E - AS) ~ ■ S(E - AS), 



(6) 



where [S(E — AS )] av = gATs{E, B) is the fully-averaged 
density of states. To calculate Q^, one can envisage 
a fully isotropic distribution of planar ^-polygons to 
obtain: 
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The second term in the numerator of Eq. ([7]) gives rise 
to a peak in Se/e at very low B, while the first term 
(present only if n w > 2) gives rise to a negative contri- 
bution to Se/e at larger B which can win over the en- 
hancement term for all values of B if Do max 3> -Domm 
{Domin, Do max corresponding to cutoffs in the distribu- 
tion of ATS energy barriers Ub)- The observations in 
Duran and BK7 indeed show a significant depression of 
e(B) for weak fields thus giving direct evidence for the 
existence of ATS with n w > 2 in the multisilicate glasses. 
Carrying out the averaging [. . .] av one gets analytical ex- 
pressions for the polarizability; the uniform average over 
orientation angles [3 must be performed numerically. 

Magnetocapacitance, Al 2 3 -BaO-Si0 2 Glass "AlBaSi-O" 
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FIG. 3: Dielectric constant's change (real part) in a magnetic 
field, measured and from the present theory, for the multi- 
component AlBaSi-0 glass at T = 118 mK (driving field 15 
kV m" 1 , driving frequency 1 kHz) [inset: data and theoreti- 
cal calculations for the BK7 and Duran glasses at 15.4 mK]. 
Data from Ref. 0|. 

Fig. 3 shows the change due to a magnetic field in the 
dielectric constant e(T, B) = e + a(T, B) for the AlBaSi- 
O glass within the present oversimplified theory, which is 
used to fit the data at 118 mK It is seen that using 
the parameters extracted from fitting this theory 12j to 
the C P (T, B) data [ll| and a value A = 20.0 x 10" 5 for 
the dimensionless prefactor A = pfP* / '(eoe r D m i n ) one 
can reproduce all the qualitative features of the puz- 
zling independence of the magnetocapacitance for this 
glass. A best fit to the 118 mK experimental data yields 
the parameters reported in Table II. The T-dependence 
of the calculated magnetocapacitance is also in qualita- 
tive agreement with the experimental data 0, |J. In a 
comparison with data for the BK7 and Duran glasses 
Q (at the much lower only available temperature of 
15.4 mK, however) the present theory reproduces the 
main features correctly. The values of Do m i n \q/e\S/\ and 
Domax\q/e\SA are very similar to those used to fit the 
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C P (T,B) data 0; the small values of D 0min /D 0max 
and of A for BK7 denote a much reduced presence of 
microcrystallites in the material. The values of D m i n are 
lower (due to the strong electric field applied) and more 
realistic than those used in the analysis of C P (T,B), in 
all cases confirming the consistency of the assumption 
\Ei\/D < 1 in this theory. 



glass 
type 


A D m in Domin\q/e\SA A)max | <?/e 1 5*A 

[mK] [K A 2 ] [K A 2 ] 


AlBaSi-0 


1.42xl(T 4 30 9.72xl0 4 1.69xl0 5 


Duran 
BK7 


0.98xl0" 5 34 1.83xl0 4 6.28xl0 5 
0.95xl0~ 5 24 4.16xl0 4 1.09xl0 6 



TABLE II: Local ATS parameters extracted from this fit 



Finally, the present simplified theory has focused 
on the weak S-field regime (up to 1 T). Given the 
large values of Do\q/e\SA thus extracted (indicating 
that small coherent clusters of some 5 to 10 NM- 
ions are involved in the magnetic-sensitive tunneling 
[13]) one can expect the low-B expression employed 
for the gap to break down around some larger field 
B* — tpo(D 0min / 'D Qmax ) I '(27r|g/ 'e\S a) above which the 
gap grows sub-linearly with 4>{B). As will be shown else- 
where, this is in turn responsible for the second enhance- 
ment of Se/e observed in the experiments for B > B* 
(B* being rather material-dependent). 

In summary, the ingredients of the present two-species 
TM together with the reasonable assumption of partial 
devitrification in the films (which can be checked through 
X-ray analysis) allows for a first good understanding of 
the puzzle of the magnetocapacitance in the cold multi- 
silicate glasses. The absence of a magnetoacoustic re- 
sponse in such glasses with ATS can be understood in 
terms of the much higher resonance frequencies of the 
NM-pockets and channels, and experiments should be 
done in such conditions. As for the PE-experiments, the 
orbital-coupling approach with the inclusion of TS-TS 
interactions has been shown to provide a partial expla- 
nation for the -B-dependence of the PE-amplitude at ul- 
tralow temperatures [l8| . The present theory, interaction 
improved, is expected to also provide an explanation for 
such data, the so-called isotope effect 0] being, possibly, 
fabrication- rather than isotope-related. Further experi- 
ments are needed: clearly, if the magnetic effects are due 
to quadrupole moments then the response should scale 
with the quadrupole-carrying nuclear concentration. If 
tunneling paramagnetic moments are involved, as is sug- 
gested by a localised TM also capable of providing a good 
explanation for the C p and e data [1 91 ] , then the magnetic 
response should scale with the Fe-concentration. In the 
present approach the response scales with the NM-ions' 
concentration and with the degree of devitrification, thus 



a larger magnetic response than thus far observed should 
be found in the best ceramic glasses, like for instance 
Ceran. 

The permission by its Authors to display Fig. 1 is 
gratefully acknowledged. 
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